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The gluon propagator in Landau gauge is calculated in quenched QCD on a large 
(32® X 64) lattice at /3 = 6.0. In order to assess finite volume and finite lattice 
spacing artefacts, we also calculate the propagator on a smaller volume for two 
different values of the lattice spacing. New structure seen in the infrared region 
survives conservative cuts to the lattice data, and serves to exclude a number of 
models that have appeared in the literature. 


1 Introduction 

The infrared behaviour of the gluon propagator is important for an und^^tand- 
ing of confinement. Several studies using Dyson-Schwinger equations tip have 
suggested that the gluon propagator diverges faster than 1/q^ in the infrared, 
and it has been claimed that this is a necessary condition for confinement. On 
the other ^nd, recent studies of coupled ghost and gluon Dyson-Schwinger 
equationsB’u indicate th^gluon propagator vanishes in the infrared, in-agree¬ 
ment with older studies HQ using other methods. Other investigations 13 have 
also suggested a massive gluon, yielding an infrared finite propagator. 

Lattice QCD should in principle be able to resolve this issue by first- 
principles, model-independent calculations. However, lattice studies have up 
to now been inconclusiveB’u since they have not been able to access sufficiently 
low momenta. The lower limit of the available momenta on the lattice is given 
by (7min = Stt/L, where L is the length of the lattice. 

Here we will report results using a lattice with a length of 3.3 fm in the 
spatial directions and 6.7 fm in the time direction. This gives us access to 
momenta as small as 400 MeV. 

2 Lattice formalism 

The gluon field can be extracted from the link variables U^{x) using 
U^{X) = 0 (^ 3 ) ^ 
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Inverting and Fourier transforming this, we obtain 




2*500 


iVM - Uli-q)) - :jTr (C/^(g) - UK-q)) 


( 2 ) 


where Uij,{q) = and A^{q) = The available momen¬ 

tum values q are given by 


q^ = 27rn^/(aL^), 


b; • - ■ 5 -bu 1 


( 3 ) 


where is the length of the box in the /r direction. The gluon propagator 
D'^tiq) is defined as 


DfM = mq)Al{-q))IV. (4) 

In the Landau gauge, the propagator is expected to have the form 

( 5 ) 

where the new momentum variable q is defined by 


2 . 5^0 

5 ^ = -sm—, 


such that at tree level, D{q^) will have the simple form 

D{q^) = l/g2 . 


( 6 ) 

( 7 ) 


3 Simulation Parameters, Finite Size Effects and Anisotropies 

We have analysed 75 configurations at /3 = 6.0, on a 32^ x 64 lattice. Using 
the value of a~^ = 1.885 GeV^ this corresponds to a physical volume of 
(3.35^ X 6.70) fm. For comparison, we have also studied an ensemble of 125 
configurations on a smaller volume of 16^ x 48, with the same lattice spacing. 
An accuracy of ^ \d^A^\^ < 10“^^ was achieved for the Landau gauge 

condition on both lattices. 

In the following, we are particularly interested in the deviation of the gluon 
propagator from the tree level form We will therefore factor out the tree 
level behaviour and plot q^D{q^) rather than D{q^) itself. Furthermore, in 
order to eliminate the most obvious source of lattice anisotropies, we will plot 
the data as a function of the momentum variable q defined in (I- 
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Figure 1: Componentwise data, for the small lattice (left) and the large lattice (right). The 
filled triangles denote momenta directed along the time axis, while the filled squares denote 
momenta directed along one of the spatial axes. 


Fig. 0 shows the gluon propagator as a function of qa for both lattices, with 
momenta in different directions plotted separately. Only an averaging over the 
three spatial directions has been performed. For low momenta, there are large 
anisotropies in the small lattice data, due to finite size effects. This can be 
seen in particular from the difference between points representing momenta 
along the time axis and those representing momenta along the spatial axes for 
qa ~ 0.4 and qa ~ 0.75. These anisotropies are absent from the data from the 
large lattice, indicating that finite size effects here are under control. 

However, at higher momenta, there are anisotropies which remain for the 
large lattice data, and which are of approximately the same magnitude for 
the two lattices. These anisotropies must be the result of finite lattice spac¬ 
ing errors in the action, which break down the continuum 0(4) symmetry to 
the hypercubic symmetry. In order to eliminate these anisotropies, we se¬ 
lect momenta lying within a cylinder of radius Ago = 2 x 27r/32 along the 
4-dimensional diagonals. The data surviving this cut are shown in fig. 

We also consider a further cut, keeping only those momenta which lie 
within a cone of 20° opening angle, directed along the diagonals and with its 
apex at the origin. This cut was found necessary to remove the low-momentum 
anisotropies on the small latticeM The data surviving this cut are shown as 
filled symbols in fig. It is interesting to see that even with this conservative 
cut, the main qualitative features of the data can still be observed. 
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Figure 2: The gluon propagator from the large lattice. Only momentum vectors with a 
distance from the diagonals of 2 units of momenta or less have been kept. The solid circles 
denote those points which also lie within a cone with an opening angle of 20°. 


4 Model Fits 


We have considered the following functional forms: 


D{q^) 

D{q^) 

where 

D{q^) 

D{q^) 

D{q^) 


(^2)1+0 _|_ ]^2 


„2 , g" + 4M2(g2) 


(92 + M2(g2))ln. 
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M(g2) =M{hi 


q2 4jvf 2 4^2 > 6/11 
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((/2 + M2)i+“ q^ + M'^ 
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(g2)i+“ + (M2)i+“ q^ + M"^ 


Z + 
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g2 + M2 


( 8 ) 

( 9 ) 


( 10 ) 

( 11 ) 

( 12 ) 


Model (| 
by Cornwall. 


is the one used by Marenzoni et ala Model (g) was proposed 
Models (0) and ( [lT| ) are constructed as generalisations of (||) 
with the correct dimension and tree level behaviour in the ultraviolet regime. 

All models were fitted to the large lattice data using the cylindrical cut. 
The lowest momentum value was excluded, as the volume dependence of this 
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Figure 3: degree of freedom for fits to model (p^. The fits are to the large lattice 

data with the cylindrical cut, with nearby points averaged. The data points are numbered 
1, ..., 142, with the most infrared point being number 1. ‘Fit Width’ denotes the number 
of points included in the fit. 


point could not be determined. In order to balance the sensitivity of the 
fit between the high- and low-momentum region, nearby data points within 
A{qa) < 0.05 were averaged. 

Model ( 0 ) accounts for the data better than any of the other models. 
Fig. ^ shows per degree of freedom for fits to this model, as a function of 
the starting point of the fit (where the points are numbered 1, 2, .. .starting 
from the most infrared) and the number of points included in the fit. The 
region of interest is the right-hand section of the plots, where the number of 
points in the fit is large and the infrared region is included. ( 0 ) gives a good 
fit (with y^/dof ~ 1) for a wide range of momenta. Only when we try to 
include the most infrared points does y^/dof go up to around 4. 

The parameter values for fits to 0 ) are shown in fig. As long as 
a reasonable number of points, including points in the infrared region, are 
included in the fit, the values vary very little when the fitting interval is varied. 
Our best estimate for the parameters is 

Z = 1.2ltlA = 1.06+1^50 = 0.7811m = O.SSla (13) 

where the errors denote the uncertainties in the last digit(s) of the parameters, 
due to fluctuations in the parameters as the fitting interval is varied. 

Fig. I shows the fit to 0) using these parameters, as well as the best fits 
to the other models considered. 
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Figure 5: Fits tofa) model (^); (b) models (dash-dotted line), d) (dotted line), (^) 
(solid line) and (|l2|) (dashed line). 
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5 Discussion and Outlook 

We have evaluated the gluon propagator on an asymmetric lattice with a large 
physical volume. By studying the anisotropies in the data, and comparing the 
data with those from a smaller lattice, we have been able to conclude that 
finite size effects are under control on this lattice. 

A clear turnover in the behaviour of q^D{q^) has been observed dX q ^ 
IGeV, indicating that the gluon propagator diverges less rapidly than \/q^ in 
the infrared, and may be infrared finite or vanishing. 

A more detailed study, including an investigation of the ten^ structure 
and a detailed analysis of different functional forms, is in progresslld The effect 
of Gribov copies is also an important issue for consideration. In the future, we 
hope to use improved actions to perform realistic simulations at larger lattice 
spacings. This would enable us to evaluate the gluon propagator on larger 
physical volumes, giving access to lower momentum values. 
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